Fuel 111 (2013) 148-156 



ELSEVIER 


Contents lists available at SciVerse ScienceDirect 

Fuel 

journal homepage: www.elsevier.com/locate/fuel 



Fast pyrolysis of coal, peat, and torrefied wood: Mass loss study 
with a drop-tube reactor, particle geometry analysis, and kinetics 
modeling 

Henrik Tolvanen *, Lauri Kokko, Risto Raiko 

Department of Chemistry and Bioengineering, Tampere University of Technology, Korkeakoulunkatu l, 33720 Tampere, Finland 



CrossMark 


HIGHLIGHTS 
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• Density and initial volume measurements are crucial in pyrolysis modeling. 
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The purpose of this study was to experimentally test and compare the fast pyrolysis behavior of torrefied 
wood, peat, and two types of coal. The experiments were conducted in a laboratory-scale drop-tube reac¬ 
tor (DTR) at a temperature range of 700-900 °C. Before pyrolysis, the sample particles were sieved with 
vibration sieves the opening of which was 100-125 pm. The initial size distribution of the sample parti¬ 
cles and their diameter evolution during pyrolysis was studied by using optical techniques. According to 
the optical measurements particle swelling during pyrolysis occurred with the tested coal types but not 
with peat or torrefied wood. With torrefied wood the particle shape changed during pyrolysis from elon¬ 
gated to spherical. The density of the samples was measured with a mercury porosimeter. The mass loss 
of the sample particles during pyrolysis was modeled with two first order models: the single-step one 
reaction kinetics model and a model in which two reactions competed to form char and volatiles from 
the virgin matter. The kinetic parameters of the reactions and the diameter evolution equation coeffi¬ 
cients were determined with both models and with all fuels. The optical measurement data from the par¬ 
ticles was used to discretize the particle size distribution. The discretized size fractions were then used in 
the model calculations instead of a mono-sized single particle approach. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Due to coal’s significant role in energy production and its harm¬ 
ful environmental impacts, new clean combustion technologies 
need to be developed if the use of coal continues. In addition to cre¬ 
ating cleaner combustion technologies, alternative fuels must also 
be found and researched to replace fossil coal as an energy source. 
Fortunately peat and torrefied wood can be considered as fuels 
capable of replacing fossil coal, to some extent even in the existing 
power plants. However, compared to coal these fuels have different 
pyrolysis and combustion properties [1,2]. The differences in the 
volatile yield and composition between these fuels can be a major 
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issue. The chemical kinetics, porosity, and particle shape behavior 
of the torrefied wood may also cause problems for example in the 
operation of a fuel burner in a boiler. Therefore, the usage of bio 
fuels in power plants and especially the modeling of these pro¬ 
cesses require more in-depth information about the fuel behavior 
during pyrolysis and combustion. 

Torrefied wood is thermally pretreated wood that has some of 
the good properties of fossil coal, such as hydrophobicity. The ben¬ 
efit of torrefied wood as a fuel is its ability to withstand outdoor 
storing without losing its heating value. Compared to untreated 
wood it is also easier to grind [3,4]. Moreover, torrefied wood 
can be considered to be carbon dioxide neutral fuel. Therefore, in 
the world of emissions trading, many energy production compa¬ 
nies have expressed more interest towards it. Nevertheless, the 
effects of torrefaction on biomass pyrolysis and combustion have 
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yet to be studied in depth. Since the torrefied wood has been 
planned as a replacement fuel for fossil coal, the kinetics studies 
should be conducted with coal references. 

Laboratory-scale testing provides useful and necessary informa¬ 
tion on solid fuel behavior during pyrolysis. This information can 
be used later on when designing larger power plants and burning 
facilities. Plenty of laboratory-scale equipment has already been 
developed for pyrolysis and combustion research. For example, a 
drop-tube reactor (DTR) can be used to simulate the temperature 
level, atmosphere, and heating rate of fluidized bed combustion 
and pulverized fuel firing; processes that are characteristic for 
power plants. Among others [5] have noticed that pyrolysis kinet¬ 
ics change substantially when the heating rate varies. Therefore, 
new advanced experimental setups should be developed in order 
to provide insight on the behavior of both fossil fuels and different 
types of biomass under conditions that resemble those of indus¬ 
trial-scale plants. As a result of the high heating rate a DTR can 
be seen as a more appropriate tool for studying solid fuel pyrolysis 
than, for example, a thermogravimeter. Experimental and model¬ 
ing studies have also been conducted in fluidized bed environment. 
Papadikis combined multiphase flows with discrete particle track¬ 
ing by incorporating an external user defined function in FLUENT 
and modeled biomass particle fast pyrolysis in a bubbling fluidized 
bed reactor [6,7]. Bruchmuller’s et al. [8] work investigated fast 
pyrolysis in a fluidized bed reactor from a particle-scale perspec¬ 
tive taking also into account mixing, segregation and entrainment 
of the biomass. However, the measurement setup used in this work 
enables the investigation of the sample particles’ actual size and 
shape before and after pyrolysis. This information is important if 
different particle geometries are taken into account in modeling. 
Among others Lu has showed in his work [9] that biomass particle 
shape has a great impact on the conversion rate. 

Several types of models have already been presented in the lit¬ 
erature and employed to describe the pyrolysis process of solid fuel 
particles [6]. However, despite the ability of these models to accu¬ 
rately describe the phenomena related to pyrolysis, strongly sim¬ 
plified kinetics models have also been used successfully [10,11]. 
If the results of these kinetic analyses are to be used in computa¬ 
tional fluid dynamics (CFD) calculations, it is clear that complex 
models are more time demanding. Therefore, simplified models 
should be favored if their accuracy is close to the complex ones. In¬ 
stead of complicated models the emphasis should be on defining 
accurate initial parameters for the modeling calculations. Thus, it 
is important to note that the initial parameters of pyrolysis model¬ 
ing, such as fuel density and diameter distribution, have a tremen¬ 
dous effect on the calculated kinetic parameters. The novelty value 
of this work lies in the complementary usage of optical techniques 
and mercury porosimeter in the sample particle shape and mass 
determination as well as in the liquid nitrogen drop-tube reactor 
collecting system. 


2. Materials and methods 

2.2. Fuel composition and sample preparation 

In this study, four different solid fuels were tested: two types of 
coal, peat, and torrefied wood. The two coal types formed a refer¬ 
ence to which the peat and torrefied wood samples could be then 
compared. The ultimate and proximate analyses of the fuels were 
conducted by Enas Co. and the results are presented in Table 1. 

It can be seen from Table 1 that there is only minor difference in 
the composition of the two coal samples. The largest fractional dif¬ 
ference was in their sulfur content. 

The torrefied wood sample material was received in the form of 
pellets. From the torrefied wood raw material 90% was softwood, 


and the rest was hardwood. The stem wood raw material was first 
chipped and then torrefied. The torrefaction was done at 250 °C in 
a continuous operated reactor. The residence time of the chips in 
the reactor was approximately 0.5 h. The mass yield of the torre¬ 
faction process was 89.6% from the original weight. The torrefied 
wood was then pelletized with a ring matrix pelletizer. 

The density of the samples used in the experiments was mea¬ 
sured with a Micromeritics Poresizer 9320 mercury porosimeter 
in the Department of Civil Engineering at Tampere University of 
Technology. The density measurements were conducted with the 
same sample particle size as used in the pyrolysis tests. The mea¬ 
sured intrinsic (skeletal) density of the sample fuel materials and 
the evaluated apparent density values of the ground sample parti¬ 
cles are presented in Table 2. The intrinsic density of the samples 
was calculated according to the mercury intrusion volume when 
all the pores greater than 6 nm were filled. The apparent density 
was evaluated with a mercury intrusion threshold of 26 pm. The 
threshold value was chosen according to the evaluated resolution 
of the already defined particle projection surface area. In the mer¬ 
cury intrusion volume curve as a function of average pore diame¬ 
ter, the 26 pm value appeared to be the point where mercury 
had filled most of the gaps between the particles and started 
intruding inside them. In this way, the evaluated particle apparent 
density and the particle volume based on image analysis could be 
used to calculate the mass of the particles. 

The coal 1 sample batches treated at 700 °C and 850 °C were 
sieved with ring sieves sized 100-125 pm and the batch at 
900 °C later on with 112-125 pm sieves. The coal 2 sample batches 
were all sieved with 100-125 pm sieves. The peat sample batches 
were sieved with 112-125 pm sieves but there were two different 
batches, one for 700 °C and 850 °C temperatures and one for 
900 °C. Finally, the torrefied wood sample batches were all sieved 
with 112-125 pm sieves. Before sieving the samples were dried in 
an oven at 105 °C for 24 h and the crushed with a Retsch ZM 200 
centrifugal mill. The sieving was conducted several times for a sin¬ 
gle sample with a vibration sieve in order to obtain a narrow size 
fraction. Between the sample preparation and the measurements, 
the sample was kept in an exicator with silica gel to remove the 
moisture from it. Therefore, in the modeling phase the particles 
were assumed to be dry when they entered the reactor. 


2.2. Drop-tube reactor 

The DTR system used in this study was designed by the first 
author. The design was based on the previous versions constructed 
at the Tampere University of Technology. The reactor consisted of 
three different modular parts: a feeding probe, a reactor part, and a 
new improved fuel collecting system. The reactor itself was an 
austenitic stainless steel tube with an inside diameter of 
26.7 mm, and with a temperature resistance up to 1300 °C. The 
reactor was covered with 2-7 Q m -1 resistance wire sets, which 
formed eight separately adjustable heating elements. The heating 
elements were insulated with an approximately 7 cm thick layer 
of kaowool. The heating zone length of the reactor was 65 cm 
and it was followed by 2.5 cm of unheated zone. The DTR and 
the auxiliary equipment are presented in Fig. 1. 

Furthermore, windows for optical measurement purposes were 
built into the reactor, and they were placed at the lower end of the 
heating zone. The center point of the windows was 53.5 cm below 
the start of the heating zone. The windows were placed as low as 
possible in the reactor to enable particle velocity measurements 
with longer residence times. As already mentioned, after the heat¬ 
ing zone end, there was a 2.5 cm section of unheated reactor tube, 
which was left between the resistance wires and the collecting sys¬ 
tem to prevent excess heating of the collecting vessel. 
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Table 1 

The results of the ultimate and proximate analyses for the tested fuels. 


Analysis 

Peat 

Coal 

1 

Coal 

2 

Torrefied 

wood 

Unit 

Ash content (815 °C) 

7.1 

13.7 

15.4 

0.2 

dry wt.% 

Volatile matter 

65.4 

34.5 

31.4 

81.9 

dry wt.% 

Sulfur 

0.14 

0.33 

1.37 

<0.02 

dry wt.% 

Carbon 

52.6 

67.8 

64.8 

53.2 

dry wt.% 

Hydrogen 

5.5 

4.6 

4.5 

6 

dry wt.% 

Nitrogen 

1.24 

2.04 

2.07 

<0.2 

dry wt.% 

Oxygen (calculated) 

33.42 

11.53 

11.86 

40.38 

dry wt.% 

Calorimetric heating 
value 

20.96 

28.05 

28.2 

21.07 

MJ/kg 

Lower heating value 

19.76 

27.06 

27.24 

19.77 

MJ/kg 


Table 2 

The measured intrinsic densities of the samples and the evaluated apparent densities 
of the sample particles. 


Sample 

Intrinsic density (kg/m 3 ) 

Apparent density (kg/m 3 ) 

Peat 

1245 

459.7 

Coal 1 

1412 

1050 

Coal 2 

1353 

1251 

Torrefied wood 

1208 

931.0 


The main function of the feeding probe was to carry the parti¬ 
cles to a wanted level inside the reactor and maintain them at a 
low enough temperature before entering the heating zone. The 
particles were inserted to the feeding probe from a feeding silo. 
Around the particle feeding tube was a water cooling jacket that 
kept the inside temperature of the probe at less than 80 °C. This en¬ 
sured that the pyrolysis process of the sample particles started 
only after they entered the reactor itself. The temperature of the 
reactor wall was constantly measured from eight points and the 
gas temperature inside the reactor was measured before every 
run. From this data, average temperature profiles were first gener¬ 
ated for the reactor wall and gas temperatures. These temperature 
profiles were then used in the energy balance equation as the gas 
and wall temperatures. The volume flow of nitrogen inside the 
reactor at 0°C was 1.5851 mol -1 , which corresponded with the 


average gas velocities of 0.1735 ms -1 , 0.200 ms -1 , and 

0.209 m s _1 at furnace temperatures of 700 °C, 850 °C and 900 °C 
respectively. The calculated Reynolds number for the reactor gas 
flow was well in the laminar zone in each of the aforementioned 
conditions. 

The collecting system was designed to quench and cool down 
the particles very rapidly. This was achieved by floating the collect¬ 
ing vessel on liquid nitrogen. The vessel was placed near the reac¬ 
tor exit and while the gases erupted out from the edges, the 
collecting vessel acted as an impactor and collected the particles. 
The evaporating nitrogen kept the atmosphere totally inert. This 
collecting system improved the accuracy of the measurements sig¬ 
nificantly. As the drop-tube reactor was constructed the mass loss 
measurement accuracy was tested with coal and biomass particles. 
A simple way of testing the system accuracy was pyrolyzing sam¬ 
ple particles in the reactor and increasing the reactor length after 
the final volatile yield was reached. If the yield stayed the same 
with greater particle falling lengths it meant that no particles were 
lost to the walls. Tests were also made in cold state to see whether 
all the particles that were fed were collected. The collecting system 
made it possible to collect all the particles that were dropped to 
the reactor, which in turn made the mass loss measurements sim¬ 
ple and accurate. 

2.3. High-speed camera 

A high-speed camera was used for taking photographs of the 
particle stream inside the reactor through the measuring windows. 
The recorded pictures were then analyzed with a computer pro¬ 
gram in order to determine the velocity of the particles. The cam¬ 
era was an AVT Marlin 145-B2 with a 1380 x 1090 resolution and 
with a black and white CCD-cell. The camera was placed in front of 
the measurement windows, and the background led-light was 
placed to the opposite side of the reactor. The falling particles 
generated a double shadow to the image because of the pulsating 
LED-light located on the opposite side of the reactor. Based on the 
distance of the shadows and the time delay between the two 
pulses, the analysis program could determine the velocity of the 
particles. However, the images had to be scaled before the 
measurements. 
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Fig. 1 . A description of the DTR on the left and the diameter measurement system on the right. 
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Before recording, the camera was focused by running some par¬ 
ticles through the reactor and adjusting the camera distance and 
aperture accordingly. Using a high-speed camera eliminated the 
need to calculate the particle residence time in the reactor with 
the particle motion equation. Instead, the velocity was first mea¬ 
sured from several stages in the reactor by moving the position 
of the particle feeding probe. Based on these measurements, veloc¬ 
ity profiles were generated for each cases and the residence time 
could then be integrated from them. An example of the velocity 
measurements and the generated fits can be seen in Fig. 2. 

Fig. 2 shows that the velocity profile of the particles was not 
even, since there was a rather large increase in the gas velocity 
after the probe tip. However, it was visually observed that this 
acceleration in the speed did not spread the particle stream. The 
acceleration in the speed resulted from the geometry inside the 
reactor. 

2.4. Determining particle diameter 

The diameter of the particles was determined in a separate 
stage after the reactor treatment. First, the particles were scattered 
on a glass plate that was luminated from below. Here, the CCD 
camera was used to take pictures of the particle projections. On 
top of the plate the particles were easier to focus than while falling 
inside the reactor. The images were analyzed with a program 
developed by PhD Markus Honkanen from Tampere University of 
Technology. The analysis program recognized the particle projec¬ 
tion outlines and calculated the pixel based projection area. In each 
diameter measurement round, altogether over thousand particles 
were detected from the images. The mass mean diameter of the 
sample particles could be determined by first calculating their vol¬ 
ume and assuming they all had the same density. Fig. 3 illustrates 
the particle projection outline recognition of the crushed and 
sieved sample particles. 

Fig. 3 shows remarkably that the vibration sieving process pro¬ 
duced mostly homogeneously sized and shaped spherical particles 
with the coal samples. Unfortunately, this was not the case for peat 
and especially for torrefied wood. The torrefied wood particles 
were highly elongated. This can most likely be explained by the fi¬ 
brous structure of wood. If the projection area of the torrefied 
wood and peat particles would have been considered as projection 
of a sphere, it would have highly exaggerated the initial particle 
mass. Therefore, the elongated particles were considered to be 
the shape of cylinders. 



Fig. 2. The measured velocity values and the generated fits for the peat sample 
particles. The zero point in the horizontal axis represents the probe tip. 


The analysis program made an ellipse fit to describe the particle 
dimensions and outlines. The dimensions of the ellipse could be 
used to estimate the sphericity of the particle. A set value of 1.5 
was used for the individual particle ellipse maximum to minimum 
length ratio to distinguish between spheres and cylinders. If the 
particle had a ratio below 1.5 when calculating its volume, it was 
assumed as sphere when calculating its volume and if equal or over 
1.5, the shape was assumed to be cylindrical. The cylinder thick¬ 
ness was calculated by dividing the projection area with the max¬ 
imum ellipse length of the particle. The maximum ellipse length 
was used instead of the minimum, since the ellipse fit seemed to 
exaggerate the minimum length for the elongated particles. With 
this approach the particle volume could be calculated more accu¬ 
rately. The presented mass mean diameters were calculated back¬ 
wards from the volume. For easier comparability, the presented 
diameters are those of a sphere with an equal volume. Therefore, 
the diameter values presented in this study describe more the vol¬ 
ume change than the actual shape change of the particles. The heat 
transfer coefficient for all particles was calculated according to the 
spherical diameter. The particle size data in turn was used to cal¬ 
culate the volume distribution of the particles. Finally, the volume 
distribution was discretized for the modeling stage. 

3. Modeling 

3.1. Modeling the mass loss chemical kinetics 

The mass loss chemical kinetics of the samples was modeled 
with two simple models: a one reaction and a two-competing reac¬ 
tions model. The one reaction model was chosen as it is commonly 
used in CFD calculations. The two competing reactions model was 
chosen to serve the need to preform solid fuel pyrolysis CFD calcu¬ 
lations with a wider particle diameter distribution. The kinetic 
parameters related to the models were calculated with an optimi¬ 
zation routine that minimized the difference between the mea¬ 
sured values and the model prediction. 

In the one reaction model there was a predetermined amount of 
volatiles in the particle that was released through one global reac¬ 
tion. The one reaction mass loss rate was expressed as follows: 

-^ = k(VMo-X DAF )m 0 (l-Ash) (1) 

The term m p is the mass of the particle (kg) at given time, m 0 is 
the initial mass of the particle (kg), t is the time (s), k is the reaction 
rate coefficient (s -1 ), andX DAF is the conversion stage (-) of the par¬ 
ticle that indicates the amount of fractional mass loss in dry ash¬ 
free (DAF) form. The term VM 0 describes the initial fractional 
amount of volatiles in the particle and it has to be determined at 
each temperature. It was also used in DAF form. The term Ash is 
the fractional amount of ash in the particle based on its initial 
dry mass. The reaction rate coefficient was of the Arrhenius form: 

k = Ae(~*&) (2) 

where A is the pre-exponential factor of the reaction rate coefficient 
(s -1 ), E a is the exponential factor of the reaction rate coefficient 
(J mol -1 ), R u is the universal gas constant (J mol -1 K -1 ), and T p is 
the particle temperature (K). In the two reaction model both reac¬ 
tions competed to form volatiles and char. Therefore, no predeter¬ 
mined amount of volatiles at a certain temperature was needed in 
the calculations. This was an advantage compared to the one reac¬ 
tion model. The mass loss rates of the two competing reactions 
were expressed as follows: 

- d Z P = ^volatiles ~ X char ~ X volatiles )m 0 ( 1 - Ash) (3) 
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Fig. 3. The particle projection outline recognition with the diameter analysis program from the sieved sample parties before pyrolysis. The green line around the particles 
illustrates the outline recognition. 


- = k char (1 - X char - X volatiles )m 0 ( 1 - Ash) (4) 

where m vo i at iies is the amount of volatiles released (kg), m char is the 
amount of virgin matter converted to char (kg), k char and k vo i a tiies de¬ 
scribe the reaction rate coefficients for the char and volatiles forma¬ 
tion reactions (s -1 ) respectively. The term X char is the fractional 
amount of DAF virgin material converted to char (-), and X V0 / atI / es 
is the fractional amount of DAF virgin material converted to vola¬ 
tiles (-). Both reactions stopped when all the virgin material had 
transformed into char or volatiles. 

3.2. Modeling the particle diameter evolution 

The diameter evolution for peat and torrefied wood particles 
during pyrolysis was modeled with the following equation: 

d = d 0 (l-X/ (5) 

where d is the diameter of the particle (m) at a given time and d 0 is 
the particle initial diameter (m). If the term p is 0, the particle diam¬ 
eter remains constant and its density decreases during pyrolysis. If 
P is negative it means that the particle size increases as the virgin 
matter conversion to volatiles advances. The maximum value of 
1/3 for p describes a situation where the density remains constant 
and the diameter decreases. This equation was suitable to describe 
the shrinkage of the peat and torrefied wood particles during pyro¬ 
lysis. Nevertheless, with the coal samples an additional term was 
introduced to the diameter evolution equation to describe the 
swelling that occurred especially at the highest temperature of 
900 °C, where also the mass loss rate was the highest. Therefore, 
the swelling term was set to be dependent on the mass loss rate 
as follows: 

d = d 0 ((l(6) 


3.3. Particle energy balance 


The following equation describes the energy balance of the 
pyrolyzing particle in the DTR: 

m p c p = 8Ah c (T gas - T p ) + A p sa(Tt all - if) + h pyr ^ (7) 


where m p is the mass of the particle (kg), c p is the specific heat 
capacity of the particle (J kg -1 1< -1 ), T p is the temperature of the par¬ 
ticle (K), t is the time (s), the term 0 describes the effect of the Stefan 
flow on the particle convective heat transfer (-), A p is the surface 
area of a spherical particle (m 2 ) with a diameter of d p , h c is the heat 
transfer coefficient (W m -2 K )* Tgas is the measured gas tempera¬ 
ture in the rector (K), s is the emissivity of the particle (-), a is the 
Stefan-Boltzmann constant (W m -2 K -4 ), T wa u is the measured wall 
temperature of the reactor (K), and h pyr is the pyrolysis enthalpy 
(J kg -1 )- The term on the left hand side describes the energy stored 
in the particle. The first term on the right hand side describes the 
heat convection to the particle from the surrounding gas, the second 
term describes the radiative heat transfer between the wall and the 
particle, and the third term describes the internal heat generation or 
consumption due to the chemical reactions taking place during 
pyrolysis. The heat transfer coefficient was calculated with the Ranz 
and Marshall [12] correlation. The heat capacity of the DAF material 
in the particle was modeled according to Merrick [13]: 


Cp : DAF 


Ru [/TeA 2 ew, 2 /Tj2 \ 2 e^Tp 

Move [Up (e T «/r P _ I ) 2 \T P ) ( e r E /T, _ if 


( 8 ) 


where M ave is the average molar mass of the elements in the DAF 
matter (kg -1 mol) and T E i and T E2 are the Einstein temperatures 
(K). The Einstein temperatures selected for the heat capacity deter¬ 
mining were 380 K and 1800 K. The heat capacity for the ash was 
calculated as follows: 


Cash — 754 + 0.586 T p 


( 9 ) 


The coefficient cp described the magnitude of the swelling. In 
both Eqs. (5) and (6) the conversion stage X was in the dry form 
in the range of 0.. .X char (l - Ash). 


The Stefan flow term 6 in the convection part of the energy bal¬ 
ance equation was used as presented in [14]. The following corre¬ 
lation applies to the Reynolds numbers up to 400 [14]: 
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Q = e (-0.6 B) 


( 10 ) 


where B is related to the total gas release from a particle (-): 


Cppi /dm p \ 

2nd p k g \ dt J 


( 11 ) 


Above c vo i is the heat capacity of the volatiles (J kg -1 K _1 ) and k g 
is the conductivity of the gas (W m _1 I< -1 ). As a result of the un¬ 
known content of the volatile gases, in this research article their 
heat capacity was assumed to be the same as that of carbon 
dioxide. 

The emissivity of the particle in the radiation part was set to be 
0.9 as proposed by [15]. Due to the small size of the particles the 
Biot number was small enough for the assumption that they were 
thermally thin. Studies suggest that biomass particle diameter of 
about 250 pm is the upper limit for a uniform temperature 
(Bi < 1) [16]. Also Papadikis et al. [7] showed in his calculations 
that with a convective heat transfer coefficient of approximately 
one third in comparison to the one in this study and with 
500 pm particles the internal temperature gradient is negligible. 
Therefore, the internal heat transfer of the particles was not taken 
into consideration in the calculations in this study. The pyrolysis 
enthalpy for all the reactions in the two models was set to be 
120 kj kg _1 (endothermic). The particle temperature was solved 
from the energy balance equation by first linearizing the radiation 
term according to [17] and then discretizing the time dependance 
of the variables. 


3.4. Discretizing the particle size distribution 

The particle size distribution was calculated for the initial sam¬ 
ples according to the analyzed diameter data. The distribution was 
discretized into ten equal-sized volume fractions. After this, the 
volume mean diameter of each fraction was calculated. The vol¬ 
ume mean diameters obtained were then used as initial diameters 
in the model calculations. Thus, the model results were volume- 
averaged. The averaging was done as follows: 

10 

F avg — ^ ^ F n Vn ( 12 ) 

n=1 

where F avg is the averaged variable (X,T p ,c p , or d p ), F n is a variable 
{XnJ P ,n, c Pt n, or d p n ) related to the nth volume fraction, and v n = 0.1 
is the volume fraction. Fig. 4 shows an example of the minimum, 
the maximum, and the volume-averaged mass loss values calcu¬ 
lated with both models. 

Fig. 4 shows that with the two competing reactions model the 
final mass loss levels of the volume fractions were not the same. 
This results from the different heating rates. With smaller particles 
the heating rate was higher and the temperature in the particles 
increased more rapidly. The higher temperature in the small parti¬ 
cles favored the reaction that formed volatiles. With the one reac¬ 
tion model, the fixed amount of volatiles resulted in even final 
volatile matter amounts between the volume fractions. Due to 
the small particle size their terminal velocity was small compared 
to the gas velocity. Thus, the velocity of the particles in the reactor 
was assumed to be the same regardless of the particle size. There¬ 
fore, there was no difference in the residence times between differ¬ 
ent volume fractions. 

3.5. Objective function and the optimization routine 

A MATLAB-based objective function was written to determine 
the kinetic parameters of the selected reactions. The objective 
function minimized the square error between the mass loss and 
diameter model predictions and the measured values. The 



Residence time [s] 

Fig. 4. The dry ash-free mass loss as a function of residence time for torrefied wood. 
The dashed lines are the maximum and minimum diameter groups and the middle 
one is the average value of all 10 volume fractions. 


optimization routine used MATLAB’s fminsearch function. The er¬ 
rors from the mass loss and diameter comparisons were summed 
in such a way that the mass loss error was emphasized ten fold 
compared to the diameter error. This was done because of the 
greater variance in the diameter measurements and because the 
mass loss was of greater importance when determining the chem¬ 
ical kinetic parameters. A lower limit of 20 kj mol _1 was set for the 
exponential factors of the reaction rate coefficients. 


4. Results and discussion 

4.1. Experimental and modeling results 

The behavior of peat and torrefied wood during pyrolysis dif¬ 
fered greatly from that of the two types of coal both in terms of 
diameter evolution and the rate of mass loss. The absolute reaction 
rates (kg s _1 ) of the torrefied wood were much higher than those of 
the coals even though their densities were similar. It was also 
noted that during pyrolysis swelling of the particles occurred with 
the coals but not with peat or torrefied wood. The volume fractions 
of the samples that were determined based on the images from the 
initial samples are presented in Fig. 5. 

The volumes of the particles for Fig. 5 were calculated according 
to the procedure described in the Section 2.4. The coal 1 histogram 
in Fig. 5 shows the effect of switching the lower sieve from 100 pm 
to 112 pm. For both coal types the diameter distribution is reason¬ 
able in relation to the sieve sizes. The peat and torrefied wood his¬ 
tograms, on the other hand show higher volume concentrations on 
significantly larger diameters than the sieve sizes. It can be ex¬ 
plained with elongated particles and the spherical equivalent 
diameter calculation. This shows that if the initial mass of the par¬ 
ticles is calculated purely based on the sieve opening sizes and an 
assumption of particle sphericity, the possibility of a major error 
increases. 

The measured mass loss in DAF form and the diameter values 
along with the mass loss fit from the two-competing reactions 
model and the diameter evolution equation fit are presented in 
Figs. 6-9. The model values were calculated up to 40 cm reactor 
length. The one reaction model results are not presented in the fig¬ 
ures. They were similar with the competing reactions model. How¬ 
ever, the error function value was always greater in all the one 
reaction model results. 
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Fig. 5. The initial sample volume fractions based on the image analysis. The temperatures in the coal 1 and peat sample legends refer to the reactor temperatures at which the 
batches were used. 


For easier comparability the measured diameter values in Figs. 7 
and 9 were calculated based on particle volumes and presented as 
diameters of spheres with equal volumes. Fig. 6 shows that the com- 
peting-reactions model was able to describe the measured mass loss 
fairly accurately with coal 1. With peat the model did not meet the 
measured values as well at the two highest temperatures quite so 
well. In this study, torrefied wood had the highest mass loss rate 
during pyrolysis. With Peat, at 850 °C and 900 °C all the measured 
values were already close to the final volatile matter. It also seemed 
that no distinct isoconversion level was reached with peat, and the 
measured mass loss values slowly increased towards the end of the 
measurements. This can be explained, for example by slower reac¬ 
tions taking place or merely by a slight error in the mass loss mea¬ 
surements. These phenomena most likely reduced the accuracy of 
the determined kinetic parameters for the peat sample. However, 


in Fig. 8 the measured mass loss values matched the model predic¬ 
tion as closely as with coal 1. 

The char yield for torrefied wood in Fig. 8 was approximately 5% 
in DAF form at 900 °C. A very similar value for woody biomass par¬ 
ticles sized 0.35 mm char yield in a DTR was reported in [18]. The 
conversion profile of the torrefied wood particles is very similar to 
that reported by Umeki et al. [19] for same sized lignocellulosic 
biomass particles in an entrained flow reactor. The residence time 
differences in reaching the maximum conversion level might be 
due to different temperature profiles in the reactors. Lehto [ 11 ] re¬ 
ported a DAF volatile yield of approximately 88% for Finnish milled 
peat tested in a DTR at 800 °C with a same particle size fraction as 
used in this work. The volatile matter amount of peat differs from 
the one measured in this work most likely due to the different 
decomposition stages of the peat samples. 
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Fig. 6. The dry ash-free mass loss as a function of residence time for peat and coal 1. 
The dots represent the measured values and the continuous and the dashed lines 
the model. 


Fig. 7. The particle diameter evolution as a function of dry mass loss for peat and 
coal 1. The dots represent the measured values and the continuous and the dashed 
lines the model. 
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Fig. 8. The dry ash-free mass loss as a function of residence time for torrefied wood 
and coal 2. The dots represent the measured values, and the continuous and the 
dashed lines the model. 
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Fig. 9. The particle diameter evolution as a function of dry mass loss for torrefied 
wood and coal 2. The dots represent the measured values, and the continuous and 
the dashed lines the model. 


The peat and coal 1 sample batches treated at 900 °C were dif¬ 
ferent than those treated at 700 °C and 850 °C. Therefore, the initial 
volume mean diameter of the particles was somewhat different be¬ 
tween them. This shows clearly in Fig. 7. Figs. 7 and 9 show that 


Table 3 

The one reaction model kinetic parameter and the diameter evolution equation 
coefficients. 


Sample 

(1/s) 

(kj/mol) 

(-) 

(-) 

Peat 

577,800 

87.46 

0.2806 

- 

Coal 1 

1230 

44.78 

-0.2019 

0.08049 

Coal 2 

1352 

42.42 

-0.1530 

0.04265 

Torrefied wood 

19,250 

56.78 

0.1989 

- 


Table 4 

The final volatile matter amount for the samples at each temperature measured with 
the DTR. 

Sample 

At 700 °C 

At 850 °C 

At 900 °C 

Peat 

0.6860 

0.7249 

0.7686 

Coal 1 

0.3459 

0.4393 

0.4656 

Coal 2 

0.3750 

0.4209 

0.4343 

Torrefied wood 

0.9067 

0.9358 

0.9489 


Table 5 

The two-competing-reactions kinetic parameters and the diameter evolution equa¬ 
tion coefficients. 

Sample 

(1/s) 

(kj/mol) 

(1/s) 

(kj/mol) 

(-) 

(-) 

Peat 

942,600 

102.7 

17,010,000 

116.7 

0.2930 

_ 

Coal 1 

276.2 

37.20 

5260 

62.97 

-0.2031 

0.07606 

Coal 2 

420.0 

38.23 

1472 

48.92 

-0.1460 

0.04323 

Torrefied wood 

10.95 

20.00 

34,550 

61.44 

0.2170 

- 


both coals swoll during the pyrolysis. However, the swelling did 
not occur with peat or torrefied wood. 

Even though 112-125 pm sieves were used the sample particle 
initial diameters were much higher in each case as can be seen 
from Figs. 7 and 9. The deviation is most likely caused by the elon¬ 
gation of the particles, especially with the peat and torrefied wood 
samples. As the particles were scattered on top of the glass after 
pyrolysis they were on their side and thus created the largest pos¬ 
sible projection from themselves. Even with the rather spherical 
coal particles this resulted in measured initial volume mean diam¬ 
eters that were slightly higher than the sieve sizes. The large initial 
diameter of the particles could also have been a result of too small 
resolution when imaging the particles. However, when the camera 
was zoomed closer, the results were more or less the same. The 
main goal in the imaging stage was to have as many particles in 
the picture with adequate resolution as possible. 

The images that were taken for the diameter analysis also re¬ 
vealed the changes in the sample particle shapes during pyrolysis. 
First of all, the coal particles seemed to remain spherical. The peat 
sample particles also seemed to maintain their elongated shape. 
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Fig. 10. The shape of the peat and torrefied wood sample particles at the end of pyrolysis. Both samples were treated at 850 °C for 40 cm in the reactor. Peat is on the left and 
torrefied wood on the right. 
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The torrefied wood sample particles, on the other hand looked like 
droplets at the end of pyrolysis. This can be seen in Fig. 10. 

The same observation could be made from the ellipse fit maxi¬ 
mum and minimum length based aspect ratios. With peat the as¬ 
pect ratio remained close to the initial value, but the ratio for 
torrefied wood decreased so that in the end it was close to one. 
The reason for this behavior could be the differences in the chem¬ 
ical composition of the peat and torrefied wood samples. For in¬ 
stance, some of the components in torrefied wood may have 
liquid intermediates. Biagini et al. [20] reported slight decrease 
in the aspect ratio of woody biomass particles. The particles in Bia- 
gini’s study were sized >0.3 mm. However, the change in the aspect 
ratio was not as high as observed in this study. Peat has a higher 
content of ligning compared to torrefied wood, which in turn has 
a higher amount of cellulose and hemicellulose in it. The high 
amount of lignin in peat could result as more durable shape struc¬ 
ture of small particles during pyrolysis. 

4.2. Fitted parameters 

The chemical kinetic parameters and the diameter evolution 
equation coefficients were calculated with the MATLAB-based 
optimization algorithm by comparing the model prediction to the 
measurements conducted. The parameters related to the one reac¬ 
tion model are presented in Table 3. 

The final volatile matter amounts for the sample fuels at all the 
measured pyrolysis temperatures are presented in Table 4. 

The volatile matter amounts in Table 4 are in dry form. They 
were determined from the measured mass loss values by calculat¬ 
ing the average of the mass loss after it had reached a clear plateau. 
The chemical kinetic parameters and the diameter evolution equa¬ 
tion coefficients related to the two-competing-reactions model are 
presented Table 5. 

The diameter evolution equation exponent /3 was negative and 
the swelling magnitude term cp was above zero for the coal sam¬ 
ples in both Tables 3 and 5. This indicates that the diameter of 
the coal particles increased due to the swelling. 

5. Conclusion 

The constructed DTR measurement system proved to be a fast 
and accurate tool for conducting mass loss and diameter evolution 
measurements with coal and biomass particles. Especially the new 
liquid nitrogen quench collecting system and the particle falling 
velocity measurement with the high-speed camera eliminated a 
large number of uncertainties from the mass loss measurements 
and from determining the particle residence times. The measure¬ 
ments indicated that the from all the four tested fuels, torrefied 
wood had the greatest volatile matter amount. With all four fuels 
the final volatile yield in the DTR was larger than the proximate 
analysis had indicated. This shows again how important it is to 
use a suitable heating rate in chemical kinetic studies in relation 
to industrial applications. The optical study of the sample particles 
demonstrated that the sample material had a clear influence on the 
shape of the particles. With fibrous material the outcome of the 
preprocessing is most likely to be elongated or cylindrical particles. 
Therefore, if the initial mass of the particles is calculated purely 
based on the sieve sizes and on the assumption of spherical particle 
geometry, there is a possibility of great error. The image analysis 


also showed that with the torrefied wood sample the shape of 
the particles seemed to change from elongated to spherical during 
pyrolysis. The volume calculation based on particle projection cou¬ 
pled with the mercury-porosimeter-based density measurement 
offered means to improve the accuracy of the particle initial mass 
calculation. The two chosen pyrolysis models were able to describe 
the mass loss behavior of the fuels fairly well. Furthermore, the 
diameter evolution equations provided means to predict the vol¬ 
ume change of the sample particles. 
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